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We study charmoinum correlators and spectral functions at finite temperature within the quenched approxima- 
tion using isotropic lattices with lattice spacing a -1 =4.86 GeV and 9.72 GeV. Although we observe some medium 
modifications of the ground state charmonium spectral function above deconfinement, we find that ground state 
charmonia (J/ip and rj c ) exist in the deconfined phase at least up to temperatures as high as 1.5T C . P-wave 
charmonia (xc) on the other hand are dissociated already at 1.12T C . 



1. Introduction 

The study of charmonium at finite temperature 
has received considerable interest ever since the 
seminal paper by Matusi and Satz pQ . The behav- 
ior of charmoinum at finite temperature was stud- 
ied using potential models with some screened po- 
tential |1I2| and it was found that excited states 
like ip' and \c ar e dissociated at T c while the 
ground state charmonia J/ip and rj c were found to 
be dissociated at 1.1T C . In last year's conference 
(Lattice 2002) lattice calculation of the charmo- 
nium spectral functions were reported |3I4| which 
indicated that in contrast to potential model ex- 
pectations ground state charmonia may exist up 
to temperatures as high as 1.5T C while \c states 
are dissociated 3 . In this paper we are going to 
present our results on the charmonium system at 
finite temperature on isotropic lattice with lat- 
tice spacing a -1 = 4.86 and 9.72 GeV. The spec- 
tral functions were extracted using the Maximum 
Entropy Method (MEM) and confirm earlier 
results on survival of the ground state and disso- 
lution of \c states. The survival of the ground 
state charmonia has been confirmed also very re- 
cently in [jy. For other studies of meson spectral 
functions at finite temperature using MEM see 



2. Numerical results 

We consider the temporal correlators of meson 
currents 



g(t,t) = J2(J&t)J(o,o)) 



(1) 



with J = q(x,T)Tq(x,r) and T = 1,75,7^,757^ 
for scalar, pseudoscalar, vector and axial vector 
channels. The low-lying states in these channels 
correspond to Xco, Vc> J/i 1 an d Xci mesons, re- 
spectively. The spectral function which we are 
interested in is related to the meson correlator by 
the integral relation 



G(t,T) 



d,Lua(uj, T) 



coshw(r- 1/(2T)) 



(2) 



*Goldhabcr Fellow 



sinhw/(2r) 

Although this relation is valid a priori only in the 
continuum, it has been shown recently that also 
in the limit of non-interacting theory it holds on 
the lattice 9 . In our simulations we used non- 
perturbatively improved Wilson fermions |10| and 
the standard Wilson action for gauge fields. The 
simulation parameters are summarized in Table 
1. The k values in this table correspond to a 
vector meson mass close to the actual J '/ip mass, 
i.e. 3.13(1) GeV at j3 = 6.64 and 3.83(2) GeV 
at (3 = 7.192 (these were estimated from the spa- 
tial correlators). In Fig. 1 we show charmonium 
spectral functions below the deconfinement tem- 
perature, T c , reconstructed using MEM for two 
different lattice spacings. One can clearly identify 
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Table 1 

Simulation parameters. T c denotes the decon- 
fincment temperature. 



(3 


a- 1 [GeV] 


K 


Size 


T/T c 


6.640 


4.86 


0.1290 


48 3 x 24 


0.75 








48 3 x 16 


1.12 








48 3 x 12 


1.50 


7.192 


9.72 


0.13114 


40 3 x 40 


0.90 








64 3 x 24 


1.50 








48 3 x 12 


3.00 



the ground state peak corresponding to J/"ip, rj c 
and Xc states. The other structures in the spec- 
tral functions are lattice artifacts as their posi- 
tion strongly depends on the lattice spacing (they 
scale roughly as 1/a). In our analysis we have 
used free lattice spectral functions from [5] as de- 
fault model (for other forms, see Ref. [5]) in MEM 
analysis. 
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Figure 1. Spectral functions at 0.75T C and a 1 = 
4.86 GeV (top) and 0.9T C and a" 1 = 9.72 GeV 
(bottom). 



Below the deconfinement temperature the tem- 
poral extent of the lattice is about 1 fm or 
larger and spectral functions can be reliably re- 
constructed. As we increase the temperature the 
temporal extent becomes smaller and smaller and 
the reconstruction of the spectral function be- 
comes more and more difficult. However, even 
without reconstructing the spectral function cer- 
tain statements about their temperature depen- 
dence can be made just by inspecting the tem- 
perature dependence of the corresponding me- 
son correlators. The temperature dependence of 
the meson correlators comes from the tempera- 
ture dependence of the integration kernel in Eq. 
121 and from the temperature dependence of the 
spectral function. In order to separate the for- 
mer effect from the temperature dependence aris- 
ing from modifications of the spectral function we 
consider the ratios G(t, T)/G recon (T) for different 
channels with G recon — J d,Lua(u>, T*) cosh(cj(T — 
l/2T))/sinh(w/2T) and T* being the smallest 
temperature available at that lattice spacing, i.e. 
T* = 0.75T C for the data set with a" 1 = 4.86 GeV 
and T* = 0.9T C for data with a" 1 = 9.72 GeV. 
Our numerical results for G/G recon are shown in 
Fig. 2. As one can see from the figure in the 
pseudoscalar channel almost no modification oc- 
curs at the level of the correlation functions up 
to temperatures 1.5T C but substantial modifica- 
tions of the correlator are seen at 3T C . In the vec- 
tor channel we see small gradual modifications of 
the correlators with increasing temperature which 
become substantial at 3T C . In contrast to this 
a strong temperature dependence is seen in the 
scalar and axial vector channels as we cross T c . 

Now let us discuss the spectral functions above 
T c . Since the physical extent of the imaginary 
time direction is rather small at high temperature 
a reliable reconstruction of the spectral function 
is difficult. In particular one finds that the re- 
constructed spectral functions are quite sensitive 
to the default model. Thus a correct choice of 
the default model is important. As we do not 
expect medium effects at energies considerably 
larger than the open charm threshold we will use 
the high energy part of the spectral function cal- 
culated below T c as a default model. At low en- 
ergy this is matched to a default model of the 
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Figure 2. G/G recon for different temperatures for 
vector and pseudoscalar channels (top) and in the 
scalar and axial vector channels (bottom). The 
data for the pseudoscalar and axial vector chan- 
nels have been shifted by a constant for better 
visualization. 

form m(u>) = miuj 2 , with mi chosen such that 
m{u>) is a continuous function in the entire energy 
interval. The vector and scalar spectral functions 
above T c are shown in Fig. 3. For comparison we 
also plot there the spectral function at 0.9T C . As 
one can see the ground state peak in the vector 
channel (the J/tjj state) exists in the deconfmcd 
phase up to the temperature T — 1.5T C while 
there is no such peak in the scalar channel even 
at 1.1T C , i.e. the Xco state is dissolved already at 
that temperature. The calculation of the errors 
shown in Fig. 3 is explained in [S]. 
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Figure 3. Spectral functions above and below de- 
confinement in the vector at (3 = 7.192 (top) and 
scalar channel at j3 = 6.64 (bottom). 
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